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Abstract 

^ . We calculate the tension a of the interface between the confined and deconfined phases 

by the histogram method in SU(3) lattice gauge theory for temporal extents of 4 
and 6 using the recent high-statistics data by QCDPAX collaboration. The results are 
a/T c 3 = 0.0292(22) and 0.0218(33) for N t = 4 and 6, respectively. The ratio (j/T c 3 shows 
a scaling violation similar to that already observed for the latent heat L. However, we 
find that the physically interesting dimensionless combinations (a 3 /L 2 T) 1//2 and aT/L 
scale within the statistical errors. 
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1 Introduction 



At the transition temperature of a first-order phase transition, a mixed state can exist 
where two different bulk phases are separated by an interface. The free energy densities 
of the two bulk phases are equal and the free energy of the mixed state is higher than 
either of the pure phases by the amount of F s = a A, where A is the area of the interface 
and a is the interface tension. If the high temperature phase transition of QCD is of 
first order, the interface tension between hadronic matter and quark-gluon plasma is an 
important parameter for the formation of quark-gluon plasma in heavy ion collisions 
and for the nucleation of hadronic matter in the early Universe; in the latter case the 
inhomogeneities generated during the phase transition could affect the light element 
nucleosynthesis ||T|. 

In recent years many methods have been introduced to compute numerically the 
interface tension in lattice SU(3) gauge theory. Because the measurement of the in- 
terface tension generally requires an extensive computing effort, so far only lattices 
with small temporal extents (values of N t ) have been used. Most of the work has been 
concentrated on the N t = 2 system |fj, |3|, |[ , some reaching N t = 4 |5], [| . 

In this work we measure the interface tension in SU(3) gauge theory with the his- 
togram method introduced by K. Binder 0. It requires high-statistics histograms of 
an order parameter at or near the phase transition point. We analyze the Polyakov 
loop histograms obtained by QCDPAX collaboration || with the recent high-statistics 
simulations on N t = 4 and N t = 6 lattices. 

Dynamics of the first order phase transition is largely governed by the latent heat 
L and the interface tension a - the "driving" and "braking" forces of the transition. 
It was already observed by QCDPAX collaboration that the dimensionless quantity 
L/T c 4 does not scale, when N t is increased from 4 to 6. As reported below, we observe 
similar behavior for the interface tension a. This makes it very difficult to estimate the 
continuum limit values of both L and a. However, when one studies physical processes 
occurring during the phase transition, a and L appear in certain combinations. We 
find that the scaling properties of these combinations are much improved. 

Let us briefly consider, as an example, the initial stage of the first order phase 
transition in the early Universe. Let us assume that the system initially exists in the 
high temperature phase and is cooled below T c with a (constant) rate C = —dT/dt. 
According to the classical nucleation theory ||, near T c the system spontaneously 
nucleates small bubbles of the low temperature phase. The nucleation rate is F dtdV oc 
T A e~ FB ' T dtdV, where F B = —LT V B + aA B is the free energy of the bubble of volume 
V B and area Ab, and T = (T c — T)/T c . When T < T c , Fb has a maximum when the 
radius is R c = 2a/ (LT); a bubble with this radius is called a critical bubble. The free 
energy of the critical bubble is given by F c = F B (R C ) = ot 2 T~ 2 T c , where 

a 2 = 16trj 3 /3L 2 T c . (1) 

If a bubble larger than R c is nucleated, it starts to expand. Expanding bubbles emit 
shock waves which reheat the supercooled matter back close to the transition temper- 
ature, thus hindering the formation of the new bubbles close to the old ones. The 
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shock waves propagate approximately with the sound velocity v s . New bubbles can 
be nucleated only in the fraction of the volume unaffected by the shock waves. This 
volume fraction is [p]] 

/(*) = i- fdt'UvKt-trnnt')]. (2) 

Jt c o 

The bubble nucleation is finished at the time t pt when f{t pt ) ~ 0. After some algebra, 
the amount of supercooling and the average distance between the nucleated bubbles 
turn out to be, 

T c -T pt ^a X - 1/2 T c , (3) 

and 

d^v s ^e x/i a X - V2 T-\ (4) 

respectively, where \ = 4 log (T c 2 /C). For the QCD phase transition in the early 
Universe [ttj, X ~ 173, T c « 150 MeV and v 2 ~ 1/3. The equations above give 
T c — T pt « 0.076a T c and d ~ 3am. The interface tension and the latent heat ap- 
pear only in the combination a in this case. 

In section |2|, we describe the histogram method: we briefly introduce the histogram 
method and apply it to the SU(3) gauge theory. Then we consider finite size corrections 
to get the interface tension in the infinite volume limit. The results are given in section 
|H and the conclusions are given in section f|. 



2 The Histogram Method 
2.1 The Probability Distribution 

In this section we introduce the histogram method |?J and apply it to the deconfinement 
phase transition of SU(3) gauge theory. Let us consider SU(3) gauge theory on a lattice 
of the size V x N t a, where V = N x x N y x N z a 3 , with a being the lattice spacing. On 
an Euclidean lattice, the temperature is given by T = l/(N t a). Close to the transition 
coupling p c (/? = 6/g 2 ), the probability distribution P(O) of a variable Q, that has 
a discontinuity at the thermodynamical limit, develops a double-peak structure. For 
example, the order parameter Q can be the absolute value of the Polyakov loop. The 
peaks correspond to the pure phase configurations and between the peaks the dominant 
contributions come from mixed state configurations where the phases are separated 
by interfaces. When the volume is increased, the peaks become more pronounced. 
The suppression of the mixed configurations is caused by the extra free energy of 
the interfaces. With the labeling of the directions such that N x and N y < N z , the 
interfaces form preferentially in the (x-y) -plane. For definiteness, in what follows we 
assume that the interfaces are oriented along the (x-y) -plane; the effects of the non- 
favored orientations are discussed in the section [2.2| . Because of the periodic boundary 
conditions, there are two interfaces with a total area 2 A = 2 N x N y a 2 . 



2 



We assume that, in large enough volumes, the distributions around the peaks are 
described by Gaussians: 

Ci exp[-(fi-a) 2 K], (5) 

where i — 1, 2 labels the two peaks, and di oc V -1 / 2 , q oc V x l 2 . Here fli and fl 2 are the 
values of an extensive variable in the confined and the deconfined phases, respectively. 
We assume that, for the time being, Q is a real- valued variable. The complication 
which arises in the case of the complex Polyakov loop will be discussed later. The 
probabilities for the system to reside in the confined or deconfined phases are given 
by, respectively, exp(— fiV/T) and exp(— f 2 V/T), except for a common normalization 
factor. Here j\ and f 2 are the free energy densities of the confined and deconfined 
phases, respectively, at the given (3 and N t : the free energy densities of both phases 
can be defined by a partial partition function summed over a subset of configurations 
which belong to each Gaussian distribution. 

We assume that in the mixed state between the peaks the two phases occupy two 
parts of the volumes, V\ and V 2 , separated by two interfaces. The probability for the 
mixed phase is 

P m (0) = c m exp[-(A V 1 + f 2 V 2 )/T — a 2A/T\. (6) 

where a is the interface tension and V\ + V 2 — V. P m depends on Q only through the 
relation 

Q = (fl^ + n 2 v 2 )/v . (7) 

Then, the whole probability distribution is given by 

P(O) = P 1 (Q) + P 2 (Q) + P m (0), (8) 

with 

Pi(fi) = a expi-fiV/T) exp[-(fi - ^f/df] (i = 1, 2) (9) 

Equation (Q) is our basic assumption on the distribution of the variable. The coefficient 
c's in equations (H) and (||) depend on the power of volume V. 

At the critical coupling (3 C in the infinite volume limit, the free energy densities f\ 
and f 2 are equal. On a finite lattice, however, in general f\ and f 2 are not equal. When 
fi and f 2 are equal, the distance between the two interfaces can vary without changing 
the total free energy: P m is a constant independent of Q. Then P(fl) becomes a sum 
of two Gaussians and a constant. On the other hand, when f\ and f 2 are not equal, 
there is no flat part connecting the two peaks. 

Let us denote the two maxima of P(fi) by p max ,i and p maXi 2 and the minimum between 
the peaks by p m i n . Let the minimum point given by 

= " < 10 > 

equal to 

tt = 7iOi+7 2 2 (11) 
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with 71 + 72 = 1, which implies V\ = 71V, V% = 72V. Then, the leading volume 
dependence, that is, the Vi dependence in the exponent, cancels in the quantity 



N 2 r> ■ 

^ _ ly t j Fmm (12) 

^ ~ 2N x N y ^ (Pmax.l) 71 (Pmax, 2 ) 72 ' 

which converges to the interface tension divided by T c 3 when V —>■ oo: 

a = alTl = lim a v . (13) 

Note that when p max ,i = Pmax,2 = Pmax, the denominator in the logarithm (p max ,i) 71 (Pmax^ 
reduces to p max - The effect of the volume dependence of the coefficient c's will be taken 
into account, when we consider finite size corrections in the next section. 

Eq. (p~3| ) gives the interface tension for each N t . To get the continuum value we still 
have to take the limit Nt — > oo (a — > 0) while keeping T = T c . 

When one applies the histogram method, one has to do some decisions which, in 
principle, are irrelevant in the infinite volume limit, but can cause systematic differences 
when working at finite volumes. Let us first elaborate further on the choice of the 
variable we adopt in the analysis. Previous studies of SU(3) gauge theory on N t = 2 
lattices have used the action density [§J or the Polyakov loop || H ; we chose the latter 
because only the Polyakov loop gives a clear separation of phases for the all cases we 
have investigated. 

The Polyakov loop Q is a complex-valued observable, and near (3 C the distribution 
p(Q) develops 4 peaks, corresponding to the confined phase at the origin and 3 degen- 
erate deconfined phases to the directions exp(z27m/3), n — 0, 1, 2. While, in principle, 
one could apply eq. (|T2"D to a complex-valued order parameter, this would require ex- 
cessive statistical quality of the data. We projected Q to the real axis by either using 
the absolute value (O a bs) or by rotating it (fi ro t), by multiplying it with exp(i27m/3), 
to the sector — tt/3 < aigQ < tt/3 and taking the real part. 

The projection Q — > f2 a t, s or fl rot causes a deformation in the probability distribution 
near the origin. Let us assume that, in a large enough volume, the original distribution 
of the peak at the origin is given by P(Q) = Cexp(— |f2| 2 /<i 2 ), where C oc V and 
d oc V~ l l 2 . The projection changes this to 

P(fiabs) = Crf27r(fiabs/rf)exp(-nL/rf 2 ) (14) 
P(Oot) = Cd3^eif(V3~n vot /d)ex P (-n 2 vot /d 2 ) (15) 

The maxima of these distributions are V2ne~°- 5 xC dfor P(fi abs ) and 3.2354 . . .xC d for 
P(f2 rot ), where C d oc V 1 ^ 2 , in accord with eq. (|5|). On the other hand, the modifications 
in the other peaks away from the origin should be small and can be neglected. 

Now let us consider what happens for <5y in these cases. After some algebras we find 
that the <5y in this case also is given by eq. ( p"2"D with minimum point eq. (|TT|). 

The second choice concerns the value of (3 used in the analysis. Using the reweighting 
technique [|TTJ we adjusted (3 so that the peaks of the histograms for f2 rot had equal 
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height. Alternatively, one may use (3 where the ratio of the weights of the peaks is 
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1/3, corresponding to 1 disordered and 3 ordered phases; for the 2-dimensional g-state 
Potts models, this method has been shown to yield the correct infinite volume (3 C up 



to exponentially suppressed corrections [Tj . However, we did not apply this procedure 
here, because for the smallest volumes available to us, the minimum point cannot be 
clearly identified to apply eq. (jjjP. Further, as discussed above, we can expect that 
eq. (|T2| ) is stable with respect to the changes in (3. Indeed, we have checked that the 
oy's calculated at j3 rnn agree with those calculated at (3q within the statistical errors. 
The final values of (3$ used in the analysis are listed in table [p. The actual distributions 
of P/3 (fl ro t) are shown in fig. [I], normalized so that p max = 1. 

Thirdly, we have to specify the way the heights are measured from the histograms. 
Due to statistical noise one should not just use the minima and maxima of the his- 
tograms in the analysis. To suppress random fluctuations we performed a third order 
polynomial fit to the histograms near the extrema and calculated the extremum values 
from the fits. The fit ranges were chosen by consulting the shape of the histogram by 
eye; again, the results remained virtually unchanged, when different ranges were used. 



2.2 Finite Size Corrections 

In this section, we discuss various finite size corrections(FSC) to eq. flTIf). Recently 
there have been a lot of interests in the analysis of FSC to the interface tension JTOl. [Tj|. 



14| ; the analysis is particularly important here because of the different geometries of the 
available lattices. The starting assumptions of the analysis are that the interfaces are 
infinitely thin and that they do not interact with each other. We take into account the 
following three possible finite size effects: Gaussian fluctuations of the order parameter 
in the bulk phases, Gaussian capillary wave fluctuations of the interfaces and the zero 
mode of the translation of interfaces to the direction perpendicular to the interfaces. 
Altogether, the ansatz becomes 

^ = E Z hulk Z zcro Z 2 cw e- 2N > N ^ N ? , i±j±k (16) 

Pmax i = i 

where i labels the direction perpendicular to the interfaces (of an area Nj N^a 2 ). Here 
Pmax stands for (p maXi i) 71 (Pmax^) 72 - in the ansatz (|16|), a is assumed to be independent 
of N x , N y , and N z , but does in general depend on N t (if we are not in the scaling 
region) . 

The bulk fluctuation term arises from the Gaussian fluctuations of the order param- 
eter. The width of the Gaussian is proportional to V _1//2 , so that the normalization is 
oc V 1 / 2 . Because p max is in the denominator in eq. (|16D , 

l/y/N x N y N z . (17) 

The zero-mode contribution which arises from the translations of the two in terface s 
to direction z with keeping the relative distance I fixed is proportional to N Z (JN X N,, 
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and the contribution from the infrared divergent part of the capillary wave fluctuation 
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is proportional to (1/ ^ N x N y ) . Together they give a factor 

Zzero OC N% , (18) 

an extra power of N z originating from the measure: dQ = (Oj — Q2)N~ 1 d£. 

The capillary wave fluctuation term excluding the infrared divergent term given by 



Bunk [131 and Caselle et al. 14 is: 



Z cw (N x ,N y ) cc J^- [ri(iN y /N x )]- 2 , (19) 



'y 



where rj is the Dedekind eta function: 

oo 

77(0) = e 2m(p/24 J] (1 - e 2mn *). (20) 



n=l 



Note that Z cw depends only on the ratio N x /N y of the interface. Not so apparent is 
the fact that Z CV1 (N X , N y ) = Z cw (N y , N x ). Because there are two interfaces, the Z t 
term is squared in eq. flip]). 
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When deriving eqs. (jig) and (|jj), we have considered the interfaces in the z direction 
only, while in eq. (|TBp we consider all interfaces including those in the x and y directions. 

In our lattices the shorter spatial lengths are equal: N = N x = N y . Rewriting 
eq. flTBD, we obtain the formula for <7y defined by eq. (|TJ): 

/V"2 / o 1 1 v 

a v (N, N z , N t ) = a-j±(c+-hgN z -- \ogN + -G(N, N z , a)) , (21) 



where c is a constant independent of N and N z and 

G(N, N„ a) = log (l + 2(^^^) 2 eM-^N z ~ N)Na/N 2 ]) . (22) 



The term G appears because we have taken into account the interfaces along (x, z)- 
and (y, z) -planes and the value of it crucially depends on the geometric feature of the 
lattice. For example, for cubic volumes G = log 3 and for long cylinders G = 0; when 
the lattice is long (N z ^> N), the effect of the planar interfaces in less advantageous 
directions is exponentially suppressed. In our case, the effect of G is significant only for 
the more cubical N t = 6 lattices. Because the analytic value of c is not so far available, 
we perform a least-mean-squares fit of eq. ( pT|) to the data with parameters a and c in 
the both N t cases. 

Higher excitations, like multiple wall (4, 6, . . . walls) configurations or configurations 
containing a large spherical droplet of one phase embedded in the other phase, can 
also contribute to the probability density between the peaks. However, they have 
considerably larger interface area and are correspondingly more suppressed. Therefore, 
we neglect their effects as well as the interaction between the two interfaces. 
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N t 


N x N y N z 


Iterations 


ftrun 


A) 


<7y (firot) 


O'V (^abs) 


7l(^abs) 


4 


12 2 x 24 


910,000 


5.6915 


5.69115 


0.0282(28) 


0.0241(27) 


0.580(37) 


4 


24 2 x 36 


712,000 


5.6925 


5.69261 


0.0308(17) 


0.0300(16) 


0.556(21) 


6 


20 3 


376,000 


5.8922 


5.89049 


0.0146(34) 


0.0123(28) 


0.569(101) 


6 


24 3 


480,000 


5.89 


5.89197 


0.0158(26) 


0.0143(22) 


0.547(62) 


6 


36 2 x 48 


1,112,000 


5.8936 


5.89402 


0.0173(25) 


0.0164(26) 


0.597(42) 



Table 1: Parameters of the runs and the results for oy. /3 mn is the value used in the 
actual simulation, and /3 is the reweighted value, oy is calculated from eq. ( [P2]) , using 
both the rotated and the absolute value of the Polyakov line histograms at (3q. 71 is 
measured from the f2 a b s histogram at /3q. 



3 Results 

The parameters of the simulations are gathered in Table [l]. The simulations were 
performed with (3 = j3 run . This was reweighted to j3o in order to obtain "equal height" 
histograms for fl TOt (fig. [I]). The values of oy [eq. ([12])] for Q TOt and fi abs , as well as 
the values of 71 for f2 abs at j3 , are also presented in Table [TJ The error analysis was 
performed with the jackknife method, with the bin sizes identical to those used in 
ref. [§. 

Let us first give the results of the interface tension in the infinite volume limit 
obtained for the rotated Polyakov line: 

a _ ^/^ _/ 0.0292(22) (N t = 4) 

<Trot -<ri± e -\ o.0218(33) (N t = 6) 1 ] 

The values of the constant c in eq. (PT|) are —1.307(43) and —1.230(58) for Nt = 4 and 
6, respectively. For N t = 6, the y 2 values are quite acceptable (x 2 /dof = 0.31). This 
implies that our assumption for FSC is reasonable. 

The results obtained for the absolute value of the Polyakov line are in remarkable 
agreement with the values given in eq. (|23|): 

f 0.0295(21) {N t = 4) 
ffabs \ 0.0218(33) (N t = 6) 1 > 

(For N t = 6, the difference can be seen only in the fourth significant digit). The values 
of the constant c are -1.267(41) (N t = 4) and -1.206(53) (N t = 6), and x 2 /dof = 0.49 
for N t = 6. We can conclude that the results do not depend on the choice of f2 rot or 
fi a bs for the order parameter. For definiteness, we use only the rotated values in the 
further analysis. 

Let us look closer at the effects of the FSC and, in particular, the function G. For 
N t = 4, the linear fit to the data gives a = 0.0317(24) without FSC, 0.0301(24) with 
FSC but with G = 0, and 0.0292(24) with the full FSC ansatz - the first two values 
agree with the last value within the statistical errors. On the other hand, when N t = 6, 
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N t 


a/T 3 


Li-ioop/T 4 


Lmcrg/T 4 


Q^l-loop 


«MCRG 


2 


0.092(4) 


2.48(5) 




0.046(3) 




4 


0.0292(22) 


4.062(85) 


2.44(24) 


0.0050(6) 


0.0084(13) 


6 


0.0218(33) 


2.395(63) 


1.80(18) 


0.0055(13) 


0.0073(18) 



Table 2: The interface tension a, the latent heat L = A(e — 3P), and the scale factor 
a = [167T(j 3 /(3L 2 T)] 1 / 2 for different iV t 's. L and a are evaluated with the 1-loop and 
MCRG beta functions. 



the corresponding values are a = 0.0185(39) (without FSC), 0.0330(39) (G = 0), and 
0.0219(39) (full FSC ansatz). In this case the function G cancels most of the correction 
introduced by the other terms in eq. (PIP- 

The Nt — 4 value above is consistent with the value a = 0.027(4) measured by Brower 
et al. || using a completely different method. Recently, Grossmann and Laursen |J 
performed a histogram method analysis for Nt = 4 and 2 lattices, using the histograms 
of the N t = 4 lattices published in refs. |§, [16]]. Their result || a — 0.025(4) is consistent 
with our result. (Although the value a = 0.040(4) which originally appeared in 0, is 
somewhat inconsistent with our value, they recalculated it.) For N t = 2, they obtained 
the value a = 0.092(4). 

The numbers above which are listed in Table |2| indicate that the ratio cr/T 3 does not 
exhibit scaling behavior for N t = 2 - 6, and not even for N t = 4 - 6. This makes the 
estimation of the continuum value of a difficult. 

A similar violation of scaling was also seen in the measurement of the latent heat from 
the same data by QCDPAX collaboration ||. Because the pressure P is continuous 
across the phase transition, the latent heat L can be calculated from the discontinuity 
of (e — 3P)/T 4 . On an Euclidean lattice this is proportional to the discontinuity in the 
average plaquette < U n >'■ 

A(e - 3P)/T 4 = -72N*(3(g)/g 3 A < U a > (25) 

where (3(g) = —adg/da is the renormalization group beta function. (3(g) can be eval- 
uated with the l(2)-loop perturbative analysis, or, since the violation of asymptotic 



scaling at these /3's is well established, with non-perturbative MCRG methods |i7 
The quantity e — 3P is preferred in lattice calculations over directly measuring e or 
e + P, because the non-perturbative corrections to the latter quantities are not known. 
For N t = 2, we use the 1-loop beta function and the plaquette gap measured by Alves 



et al. pq| , and for Nt = 4, 6 we use the results of QCDPAX collaboration ||, evaluated 
with the 1-loop beta function and with the non-perturbative beta function. The values 
are collected in Table |2|. 

From L and a we can now calculate the ratio a 2 = 167rt 3 /(3L 2 T), which is also 
shown in Table |2|. We observe that when N t = 4 - 6, the scaling violations seem to be 
absent in a, at least within the accuracy reached in this work. This holds both for the 
1-loop and MCRG corrected values. On the other hand, the N t = 2 value is completely 
off the mark when compared with the higher N t results. 
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N t 


(aT/L) i_i oop 


(aT / L)mcrg 


2 


0.037(2) 




4 


0.0072(6) 


0.0120(15) 


6 


0.0091(14) 


0.0121(22) 



Table 3: The ratio aT/L, evaluated both with the 1-loop and MCRG improved beta 
functions. 



Inserting the N t = 6 value of ckmcrg 1° t ne eqs. (^) and (|j), we find that the degree 
of supercooling at the deconfmement phase transition in the early universe and the 
average distance between the nucleation centers (~ the scale of the inhomogeneities 
generated during the transition) are given by 

(T c - T pt )/T c = 5.6(1.4) x 10" 4 , (26) 

and 

d = 22(5) mm, (27) 

respectively. Another interesting quantity is the ratio aT/L, which is relevant when 
one calculates the hydrodynamical properties of the propagating phase transition front 



19|1 . These values are shown in Table y. For iV t = 4 - 6, the ratio aT/L also appears 



to be scaling within the statistical errors. 



4 Conclusions 

We have measured the interface tension a between the confining and deconfining phases 
of SU(3) lattice gauge theory by the histogram method using the data from the high- 
statistics simulations on Nt = 4 and 6 lattices by QCDPAX collaboration Mi. The main 
advantage of the histogram method is that it offers a direct way to measure numerically 
the interface tension from order parameter distributions. Hence, the results of ordinary 
finite temperature simulations near the transition point can be used. We also made a 
careful analysis of finite size corrections to the interface tension in order to estimate the 
infinite volume values. We found our ansatz for finite size corrections is quite reasonable 
in the sense that our data are well fitted to the formula with quite acceptable \ 2 values. 

The results, a/T c 3 = 0.0292(22) and 0.0218(33) for N t = 4 and 6, respectively, exhibit 
a scaling violation qualitatively similar to the latent heat L ||. Since both the latent 
heat L and the interface tension a characterize the strength of the first order transition, 
it may not be surprising that we observe similar scaling violations. Thus, it is plausible 
that when one takes a suitable combination of L and a, scaling becomes much better. 
Indeed, we observed that the physically relevant combinations a 2 oc a 3 T/L 2 and a/LT 
scale within the statistical errors. Clearly, scaling for the dimensionless combinations 
and the lack of scaling for L and a themselves are not satisfactory. However, we can 
use the combinations that do scale to predict the continuum physics. 
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Nevertheless, one has to be very cautious when applying these numbers to cosmol- 
ogy. We have neglected the effect of quarks, which may change the character of the 
transition completely. However, we believe that our results offer some intuition about 
the expected physical scales relevant for the QCD phase transition. 
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Figure captions 

Figure 1: The probability distributions of the rotated Polyakov loop, normalized to 
Pmax — 1: (a) for N t = 4 with two different volumes; (b) for N t — 6 with three different 
volumes. 
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